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Abstract 

In this paper we study the appearance of branches of relative periodic orbits in Hamil- 
tonian Hopf bifurcation processes in the presence of compact symmetry groups that do 
not generically exist in the dissipative framework. The theoretical study is illustrated with 
several examples. 



1 Introduction 

Let (V,lu) be a symplectic vector space and G be a compact Lie group acting linearly and 
symplectically on V . Let now be a one-parameter family of G-invariant Hamiltonians h\ £ 
C°°(V) G such that for each value of the parameter A, the origin is an equilibrium of the 
associated Hamiltonian vector field, that is, dh\(0) = for arbitrary A. In this paper we will 
study the nonlinear implications of the following linear behavior: suppose that there is a value 
of the parameter A Q and a pair of eigenvalues ±iv in the spectrum of the linearization at zero 



of the dynamics induced by the Hamiltonian vector field Xh x that behave as in Figure 1.1 
when we move the parameter A around A - Such a behavior in the parametrical motion of the 
eigenvalues is usually referred to as Hamiltonian Hopf bifurcation [vdM85|, denomination 



that we will use here, even though it also appears in the literature as 1 : — 1 resonance, 1 : 1 
non-semisimple resonance, and Krein collision. The reference to the Hopf bifurcation comes 
from the analogy with the codimension one non-conservative case in which a one-parameter 
family of vector fields has a pair of eigenvalues that cross the imaginary axis at a critical 
value of the parameter (the "classical" Hopf bifurcation). The case of G-equivariant vector 
fields (G compact) has led to the successful theory of Hopf bifurcation with symmetry which 



was initiated by [GoS85] and which was exposed in its most achieved form in [Fi94| (see also 



|ChL00] for a comprehensive exposition) 
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Figure 1.1: Motion of eigenvalues in a Hamiltonian Hopf bifurcation. 



The history of the Hamiltonian Hopf bifurcation in the non-symmetric setup is very long 
and we shall not attempt to survey it here. We just refer to [ McyS71| , Mcy86 , vdM85 , vdM96 , 
Bri9C , GMSD95 | and references therein for discussions. 



The only works that we know of dealing with the Hamiltonian symmetric case are | vdM9C ] 
and HKMS96B . In the first paper it is shown that one can apply the non-symmetric results 
on Hamiltonian Hopf bifurcation on some of the fixed point spaces corresponding to isotropy 
subgroups of the symmetry of the system, provided that certain dimensional restrictions are 
fulfilled. [ KMS96| | studies branches of (stable) three-tori that can be obtained out of a Hamil- 
tonian Hopf bifurcation process with a symmetry given by the semi-direct product of with 
T 2 x S 1 . See also [pri90a|| . 



Natural dynamical elements that show up in the study of systems that present a continuous 
symmetry group G are the so called relative equilibria (RE) and relative periodic orbits 
(RPOs), that is, motions that project onto equilibria and periodic orbits in the quotient space 
V/G, respectively. In our work we will see that whenever a Hopf-like motion of eigenvalues 
occurs in a Hamiltonian system with symmetry, one can prove the existence of periodic and 
relative periodic motions at the non linear level, for values of the parameter nearby A , whose 
number we will estimate at each energy level. The existence of periodic motions, in the presence 
of some dimensional restrictions that we have eliminated, was something known to the above 
quoted authors. As to the relative periodic orbits, they are found in those papers only after a 
reduction has been performed that makes the problem equivalent to that of searching periodic 
orbits in the relevant quotient space. Since this reduction cannot always be carried out in a 
straightforward manner we will follow an approach in which the existence of relative periodic 
orbits is proved in the original space V . 

Our approach to this problem will be based in the combined use of five tools: 



(i) Reduction method of Vanderbauwhede and van der Meer [ VvdM95|: it allows us to substi- 
tute the search of periodic and relative periodic orbits by the search of relative equilibria 
of a S^-symmetric associated Hamiltonian system (usually referred to as normal form). 
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(ii) Generic structure of the generalized eigenspaces corresponding to the colliding eigenval- 
ues [DMM92]: it determines the most plausible reduced space in which we should work 



after applying (i). 



(iii) Equivariant Williamson normal form | MD93|| : it is used to normalize the linear term of 



the equation that defines the S 1 ~relative equilibria that we are looking for. 

(iv) Lyapunov- Schmidt reduction of the finite dimensional equation that defines the 5 1 - 
relative equilibria and formulation of the problem in terms of a bifurcation equation 



of gradient nature [CLOR99| with very specific equivariance properties, 
(v) Solution of the bifurcation equation using either topological or analytical methods. 
The paper is structured as follows: 

• In Section |2| we briefly review the abovementioned tools, set the notation that will be 
used throughout the paper, and explain in detail the hypotheses under which we will 
work, along with their implications. The expert can skip this section and use it just as a 
glossary. 



Section H is devoted to Theorem 3.3, which provides a lower estimate on the number of 



periodic and relative periodic branches that bifurcate from the origin if there is a collision 



of eigenvalues as in Figure |L1 



In Section 4.1 we study a system of two nonlinearly coupled harmonic oscillators in the 
presence of a magnetic field, that will lead us to the consideration of the general case of 
the Hamiltonian Hopf bifurcation in the presence of a 0(2) symmetry. We will see that, in 
contrast with the dissipative case, the 0(2)-symmetry in a Hamiltonian Hopf bifurcation 
process gives rise to the appearance of numerous relative periodic motions. This example 
will also show that, in general, the topological methods utilized in Theorem |3.3| are not 
powerful enough to detect all the periodic and relative periodic elements of a particular 
system with a given symmetry, that is, the generality of this result is paid with its lack 
of sharpness. This circumstance will motivate a more hands-on approach to the problem 
in Section [4.3| , where we will see that under additional hypotheses on the group action, 
sharper general results can be formulated that give account of all the dynamical richness 



evidenced in the example in Section 4.1 



In Section |4.4| we use the previous results to show the existence of RPOs in Hamiltonian 



Hopf phenomena with spheric symmetry. 

For the sake of the clarity in the exposition, the proofs of some of the technical results 
needed in the main theorem are relegated to an appendix (Section |J) at the end of the 
paper. 
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2 Preliminaries and setup 

All throughout this paper, our discussions will mostly take place in a finite dimensional sym- 
plectic vector space (V, oS) on which the compact Lie group G acts linearly and canonically, that 
is, respecting the symplectic form uj. We will be interested in a one-parameter family of G— 
equivariant Hamiltonian vector fields X^ x , induced by the family of G-invariant Hamiltonians 
h x € C°°(V) G , A£l, such that: 

(HI) h x (0) = and dh x (0) = for all A. 

(H2) There is a value A of the parameter A for which the G-equivariant infinitesimally sym- 
plectic linear map A\ o := DyX^ x (0) is non singular and has ±i^ in its spectrum 
(fo^O). 



In the following paragraphs we introduce some tools and notations that will be used all 
throughout the paper. 

The resonance space Let (V, uj) be a symplectic vector space. It is easy to show that there 
is a bijection between linear Hamiltonian vector fields on (V,u>) and quadratic forms on V. 
Indeed, if A : V — > V is an infinitesimally symplectic linear map, that is, a linear Hamiltonian 
vector field on (V,uj), its corresponding Hamiltonian function is given by 

Qa(v) := -uj(Av,v), for any v 6 V. 

Also, since A belongs to the symplectic Lie algebra sp(V), it admits a unique Jordan 
Chevalley decomposition [Hu72, VvdM95| ] of the form A = A s + A n , where A s € sp(V) 
is semisimple (complex diagonalizable) , A n £ sp(V) is nilpotent, and L4 s ,^4 n ] = 0. If iu is one 
of the eigenvalues of A G sp(V) and T Uo := we define the resonance space U Uo of A with 
primitive period T Vo as 

U Uo := ker (e AsT "° - I) . 
The resonance space U Uo has the following properties (see |Wil36| , |GoS87| , |VvdM95| ): 

(i) U Vo is equal to the direct sum of the real generalized eigenspaces of A corresponding to 

eigenvalues of the form ±ikv , with k € N*. 

(ii) The pair (U Vo ,u}\u„ a ) is a symplectic subspace of (V,uj). 

— A i 

(iii) The mapping 9 E S i— > e v ° B \u Vo generates a symplectic S -linear action on (U Vo , ui\u Vo ), 

whose associated equivariant momentum map will be denoted by J : U Uo — > Lie(S' 1 )* ~ R. 

(iv) If (V, uj) is a symplectic representation space of the Lie group G and the Hamiltonian 
vector field A is G-equivariant (equivalently, the quadratic form Qa is G-invariant), 
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then the symplectic resonance subspace [U Vo ,u\u v ) is also G-invariant (this follows from 
the uniqueness of the Jordan-Chevalley decomposition of A, which implies that if A 
is G-equi variant, so is A s ). Moreover, the 5 1 and G actions on (U Vo , u>\u Vo ) commute, 
which therefore defines a symplectic linear action of G x S 1 on U Uo . See the Appendix 
(Section ||) for a sketch of the proof of some of these facts. 



The normal form reduction [ vdM85| , |vdM9C , VvdM95] Let (V, u),h\) be a A-parameter 



family (A € A, where A is a Banach space) of G-Hamiltonian systems such that h\ o (0) = 0, 
d^Ao(O) = 0) an d the G-equivariant infinitesimally symplectic linear map A := DXh Xo (0) is 
non singular and has ±ii> as eigenvalues. Let (U Uo ,tu\u„ ) be the resonance space of A with 
primitive period T Vo . For each k > there are a C -mapping xj) : U Uo x A -> V and a C k+1 - 
mapping h\ : U Uo xA^M such that tp(0, A) = 0, for all A £ A, Du Vo i()(Q, A ) = Iu Ua , and h\ is a 
G x 5 1 -invariant function that coincides with h\ up to order k+1. The interest of normalization 
is given by the fact that one can prove [ VvdM95| , Theorem 3.2] that if we stay close enough to 



zero in U Va and to A € A, then the 5 1 -relative equilibria of the G x S^-invariant Hamiltonian 
h\ are mapped by ip(-,\) to the set of periodic solutions of (V, u, h\) in a neighborhood of 
€ V with periods close to T Uo . Hence, in our future discussion we will substitute the problem 
of searching periodic orbits for (V,ui, h\) by that of searching the 5 1 ™relative equilibria of the 
G x 5 1 -invariant family of Hamiltonian systems (U Uo ,u\u v , h\), that will be referred to as the 
equivalent system. Note that the properties of tp imply that 

A := A\ v „ o = D v X hXo (0)\ Uuo = D v X hxo]Uvn (0) = D v X^(0). (2.1) 



Generic structure of the resonance space U Uo and canonical form of the symplectic 
pair (u\u vo ,A) After the remarks previously made, we know that the resonance space U Ua is a 
G x S^-symplectic vector space. The decomposition of U Vo into G x S 1 -irreducible subspaces 
that can be generically expected when the eigenvalues behave parametrically as in Figure 



1.1 



has been studied in [DMM92], where the authors concluded (Proposition 6.1 (3)) that the only 
possibility is 

U Vo =U l ®U 2 , (2.2) 



where U\ and XJi are complex dual irreducible subspaces of U Uo in the sense of HMRS88| , Theorem 
2.1]. In all that follows we will assume that we are in this generic situation. 



Once we know the decomposition ( |2.2| ) of U Vo into irreducibles, the equivariant version of 
the Williamson normal form of Melbourne and Dellnitz guarantees iMD93| , Table 2] that there 



is a basis of the vector space U Uo in which the simultaneous matricial expressions of uj\u and 
A are either 



(i) 



A= ( U °^ n v j j and ^| r „ : Ji„. or (2.3) 
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(ii) 

where 2n = dim U\ = dim XJi , l2n is the 2n-dimensional identity matrix, and J2n. is denned as 

hn= {l o)- 

Given that the treatment of cases (i) and (ii) is completely analogous, we will focus in all 
that follows in expression ( |2.3D , Moreover, whenever our family of G-Hamiltonian systems 
falls into the generic situation described in this paragraph, we will say that it satisfies the 
condition (H3). For clarity and future reference we state this condition explicitely: 

(H3) The resonance space U Vo corresponding to the eigenvalues ±iv a splits into two complex 
dual G x 5 1 -irreducible subspaces. This condition is generic. 

Relative equilibria and critical points Being consistent with the notation previously in- 
troduced, let A = A s + A n be the Jordan-Chevalley decomposition of A G sp G (U Uo ). We will 
denote by J : U Va — > Lie(5 1 )* ~ R the equivariant momentum map associated to the symplectic 

5 1 -linear action defined by (9,v) >— > e^^v, 6 G S 1 , v G U Uo . Also, for any £ G L\e(S 1 ) ~ R 
and any v G U Vo , we will write 3^(v) := J(u)£. The linearity of the action implies that, for 
any £ G Lie(5 x ) ~ R and any v G U Uo , the momentum map J is uniquely determined by the 
expression 

j((v) = ~u\ Uvo (£-v,v), 

where the dot in £ • v means the associated representation of the Lie algebra Lie(S' 1 ) on U Va 
through the S^-action. More specifically, 

For future reference we note that this relation implies that 

d 2 3(0)(v,w) = u\ Uua (A s v,w), for any v,weU Uo , (2.5) 
which in the basis used to write ( |2.3| ) admits the following matricial expression: 

d2j <°» = ( -I k o ) • ^ 
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A very interesting feature of the Hamiltonian framework is that the search for relative equi- 
libria reduces to the determination of the critical points of the so called augmented Hamil- 
tonian \ AM78 |. In the particular case that we are dealing with, this remark translates into 
saying that the equivalent system (U Uo ,uj\u ,h\) has a S^-relative equilibrium at v € U Uo 
(which represents a periodic orbit of the original system (V, u>, h\) with period near T Va ) if and 
only if there is an element £ G Lie(5' 1 ) for which 

d(£-J%) =0. (2.7) 



Whenever we find a pair that satisfies ( |2.7| ), we will say that v is a relative equilibrium 

with velocity £. 

Expression (2.7) can be written as a gradient equation, which will be exploited profusely 
in our subsequent discussion. Indeed, let (•,-) be a G x S^-invariant inner product on U Va 
(always available by the compactness of G x S 1 ). For any v € U Uo , we define the gradient 
{h\ — 3^){v) as the unique element in U Va , such that for w £ C/^ arbitrary 

d(/^ - J%) • w HVa^ - J% W 

Also for future reference, we recall that the linearization A\ = DyX^ifS) of at E V, is 
a linear G-equivariant Hamiltonian vector field with associated quadratic Hamiltonian function 
Qx given by 

Q x (v) = ^d 2 h x (0)(v,v), 

that is: 

i Ax uj = dQ\. 

The restriction A of A\ o to t/j, is of course also Hamiltonian but in this case, by (|2.1|), the 
associated quadratic Hamiltonian function can be expressed in terms of the Hessian at of the 
equivalent Hamiltonian h Xo associated to h Xo . Indeed, 

'^\u va =dQ Ao , (2.8) 

where, for any v G U Vo , 

Q A » = id 2 ^ o (0)(«,u) = ^d 2 /^(0)(^). (2.9) 



If we write ( |2.8[ ) using the basis that produced the canonical form fl2.3|), the equality ( J2.9| ) 
guarantees that 

d 2 /^(0) = d 2 h Xo (0)\u Xo = ~h n A = ( ° ) . (2.10) 



2/i 
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Invariant splitting of the resonance space U Va Using expressions ( |2.5| ) and ( p. 10 ) we can 



immediately construct a very convenient splitting of the resonance space U Vo : let L : XJ Vo — > U Uo 
be the linear map denned by (L(v),w) = d 2 (h\ — 3 Va )(0)(v,w), for any v,w G U Uo . Using 
expressions ( |2 . 5[ ) , ( |2.6| ), and ( |2,10p we can write, using the basis introduced in ( |2.3| ), that we 
will use in all that follows: 



L 





-hn 



Since the linear map L is G x S^-equivariant and self-adjoint we can split U Uo = Vo ® V\ as 
the direct sum of the two G x S 1 — invariant subspaces, 



U n := kerL 



aeR 2n ^, Ui:=ImL 



b G R 2n } . (2.11) 



Since by hypothesis (H3) we are in the generic situation, the resonance space U Uo splits as the 
sum of two complex dual irreducible subspaces with respect to the Gx S^-representation [ DMM92] 



Given that by construction Vq and V\ are G x S^-invariant and have the same dimension, the 
G x S^-representations on Vo and V\ are necessarily complex irreducible. We describe more 
precisely the interplay between the decomposition U Vo = Vq V\ and the G x S' 1 -action in the 
following lemma. 

Lemma 2.1 In all the matricial statements bellow we assume the use of the basis of the canon- 



ical form ( pOy . 

(a) Let g G G x S 1 arbitrary and v = v o + V\ G U Uo , with vo G Vo and v\ G V\. Then, there 
exists a orthogonal matrix A g such that [A g ,§2n] = and 



9 ■ v 



A g W vo 

Ag J V «i 



(b) The inner product on U Uo that takes the Euclidean form when expressed in the coordinates 
corresponding to the basis used to write the canonical form is G x S 1 -invariant. 



Proof (a) The G x S^-invariance of the spaces Vo and V\ implies, for any g G G x S , the 
existence of two invertible matrices A g and B g such that 



g ■ v 



A g \ ( vo 
B a \ Vl 



for any v = vo + v\ G U Uo , with vo G Vo and v\ G V\. Given that by hypothesis the family 
of Hamiltonians h\ is G-invariant, the linearization A in ( |2.3| ) of the vector field X^ x at the 
origin is necessarily G x S^-equivariant and, consequently 

A g \ f V hn hn \ = ( Vohn hn \ ( A 9 

B q )\ vj 2n ) 1 vj 2n Mo B a 



s 



which implies that A g = B g and that A g S2n = $2nA g , for any g € G x S . We now see that 
the matrices A g are orthogonal: given that the G x S^-action is canonical, we have that 



A T g \( -I 2n \ ( A g \ f -I 2n 



o ; v hn o y v o a 9 j \i 2n o 



(2.12) 



This equality guarantees that A T g A g = A g A T g = hn, as required. 

(b) Let (•, •) be the inner product in the statement, that is, for any v, w 6 U Vo whose coordinates 



in the basis of fl2.3|) are v = (vi,... ,f4 n ) and w = (wi,... ,w^ n ) we have that (v,w) 



viw\ + . . . + Vi n Wi n . Notice that, also in this basis, the inner product (•, •) can be expressed as 

(v,w) = 0j\u Vo {l4nV,w). 

Before we proceed notice that expression ( p,12| ) together with the orthogonality of A g implies 
that for any v G U Uo and any g € G we have that $4 n g ■ v = g ■ J^nV. If we put together this 
fact with the canonical character of this action we obtain the invariance of the inner product. 
Indeed, for v,w £ U Uo and g G G x S 1 arbitrary, we have that 

{g-v,g-w}= u)\ Uvo (I 4n g -v,g ■ w) = u\ Uvo (g ■ S in v,g ■ w) = u\ Vvo (J An v,w) = (v,w), 

which concludes the proof. ■ 

Remark 2.2 In all our subsequent discussions we will use the inner product presented in the 



previous lemma and the basis of the canonical form (2.3). ♦ 



The quadratic part of the Hamiltonian and a final generic hypothesis The complex 
irreducibility of the G x S^-action on Vq implies [ |GSS8§| , Lemma 3.4] that if VgxS 1 (Vb) denotes 
the ring of real G x S^-invariant polynomials on Vq, one can choose a basis {iq, ... ,-F/} of 
'PgxS 1 (Vqi Vq), that is, the finite type "Pcxs^V^-module of G x S^-equi variant polynomial 
mappings of Vq into itself, such that 

Ft = hn 

F 2 =hn (2.13) 
degF fc > 1 Vk> 2. 

Analogously, one can choose a Hilbert basis {6q, . . . , 9 r } of the module "PgxS 1 (^b)> such that 

~ W = HI 3 (2 14) 

deg^ fc > 2 Vfc > 1. 1 ' j 

In particular, the GxS^-invariance of the Hamiltonians h\ of the equivalent system (U Uo , oj\u v , h\) 
implies that for each A, the second derivative d 2 h\(0), considered as a linear map d 2 h\(0) : 
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Vo© Vi — > Vo© V\ is G x 5 1 -equivariant. At the same time, since it is a Hessian, it is symmetric 
and therefore there are functions a,p,T,ip € C°°(IR) such that: 



d 2 M0) 



a(X)l 2n T(X)I 2n + V(A)J 2 n 

r(A)I 2 „ - V(A)J 2 n p(A)I 2 n- 



where, by ( 2.10 ), we have the following initial conditions: cr(A ) = 0,p(A o ) = — l,r(A 
and ip(X ) = v . In all that follows we will assume the following generic hypothesis: 



(2.15) 
) = 0, 



(H4) The one-parameter family of G-Hamiltonian systems (V,lo, h\) satisfies that cr'(\ ) ^ 0, 
where cr(A) € C°°(M) is the smooth real function introduced in ( |2.15 ). 

The previous generic hypothesis that will be of much technical importance in what follows is 



related, as we show in the following lemma, to the motion of eigenvalues depicted in Figure 1.1 



Lemma 2.3 Let (V,u,h\) be a one-parameter family of G -Hamiltonian systems satisfying 
hypotheses (HI) through (H4). Then, there is a part of the spectrum of the linearization 
of the Hamiltonian vector fields X^ x at zero that behaves as in Figure LA as we move the 
parameter A. 

Proof We first compute the eigenvalues of A\, that is, the restriction of A\ to the symplectic 
subspace (U Ua , oj\u Vo )- This vector field is Hamiltonian. More specifically: 



dQx, 



where, for any v G U, 



Vo 1 



Q\(v) = ~d 2 h x (0)(y, v) = ^d 2 M0)( W , v). 



(2.16) 



(2.17) 



If we use in (2.17) the specific form for the symplectic form uj introduced in ( |2.3| ) we obtain 
that: 



A X = J 4 nd 2 /l A (0). 



We now use ( |2.15 ) and obtain that: 



T(\)hn + 4>(\)hn ~p(\)hn 

a(X)I 2n r(A)I 2n + ^(A)J 2n 



(2.18) 



(2.19) 



We compute the eigenvalues of this matrix using the well-known fact [Hal74, p. 102, ex. 9] 
that if A, B, C, and D commute, then 



det 



A B 
C D 



det(AD - BC). 
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Indeed, using this relation, we see that the characteristic polynomial of A\ is 

i f( A _ nJ x _ ( (-T(A) - /i)I 2n + V(A)J 2 n -p(A)I 2n \ 

[AX M 4nj " V ^(A)I 2 n (r(A) - M )I 2n + V(A)J 2 n J 

= det [(^ - r(A) 2 - ^(A) 2 + p(X)a(X))h n - 2^(A)J 2n ] 

- Hpt f ^ " r ( A ) 2 " ^( A ) 2 + ^M^) 1 " 2^(A)I n 

V -2^(A)I„ (M - r(A) 2 - V(A) 2 + p(A)a(A))I„ 

= det [(( M 2 - r(A) 2 - ^(A) 2 + p(X)a(X)) 2 + (2 / uV(A)) 2 )I„] 

= (((/J 2 - r(A) 2 - ip(X) 2 + p(A)a(A)) 2 + 4/^(A) 2 ) n . 

Consequently, the eigenvalues of A\ are 

MA) = ±V r ( A ) 2 " KAMA) - V(A) 2 ± 2|^(A)| v / p(A)a(A) (2-20) 

We are now in position to show that hypothesis (H4) implies the motion depicted in 
Figure [O]. Indeed, assume that a'(X ) / 0. Let /i(A) = p(X)a(X) - r(A) 2 . Note that 
/i(A ) = and f{(X ) = —a'(X ), that by hypothesis is different from zero. This implies 
that the function f\ changes sign at A D . More explicitely, suppose that <t'(A ) < (the case 
a'(Xo) > is completely analogous); in that case /i(A) < for A < A and /i(A) > for 
A > A D . Since we have that: 

MA) = ±\l-h(x)-i>(xy±2\^{x)Wh(^ = V-(VX(A)±IVKW. 

then, by considering the cases A < A D , A = A , and A > A in the previous expression taking 
into account the changes of sign in /i, we obtain the evolution of eigenvalues illustrated in 
Figure [□]. ■ 

Remark 2.4 The generic hypothesis (H4) is a sufficient but not necessary condition for ob- 
taining a behavior of the eigenvalues as in Figure that is, such an evolution can take place 
even for systems in which cr'(A ) = 0. ♦ 



3 Hamiltonian Hopf bifurcation and relative periodic orbits 

The main goal of this section is the statement and proof of a result that will provide and estimate 
on the number of relative periodic orbits of a one-parameter family of G-Hamiltonian systems 
(V, u>, h\) that satisfies the hypotheses (HI) through (H4), formulated in the previous section. 

We will begin by introducing some classical definitions that will make more explicit some 
of the concepts used in the previous paragraphs. 

It appears very frequently in examples dealing with symmetric families of Hamiltonian 
systems that the canonical symmetry group G contains a continuous globally Hamiltonian 
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symmetry: suppose that G contains a Lie subgroup H of positive dimension. We say that the 
canonical action of H on V is globally Hamiltonian when we can associate to it an equivariant 
momentum map K : V — > f)* which is defined by the fact that its components := (K,£) G 
C°°(R), £ G f), have as associated Hamiltonian vector fields the infinitesimal generators of the 
action 

exp i£ • v £ G f), v G V. 

t=o 

Definition 3.1 Let (V,to,h) be a Hamiltonian system with a symmetry given by the canonical 
action of the Lie group H on V . The point v (zV is called a relative periodic point (RPP), 

if there is a r > and an element g G H such that 

F t+T (v) = g ■ F t (v) for any t G R, 

where F t is the flow of the Hamiltonian vector field . The set 

7 (f) := {F t {y) | t > 0} 

is called a relative periodic orbit (RPO) through v. The constant r > is its relative 
period and the group element g G H is its phase shift. 



Proposition 3.2 Let (V,uj,h) be a Hamiltonian system with a globally Hamiltonian symmetry 
given by the canonical action of the Lie group H on V with associated momentum map K : 
V — > f)*. If the Hamiltonian vector field X h _^, £ G f), has a periodic point v G V with period 
t, then the point v is a RPP of Xh with relative period r and phase shift expr£. 

Proof Let Ft be the flow of the Hamiltonian vector field Xh and Kt(v) = expi£ • v that of 
Xj£$. By Noether's Theorem: 

[Xh,XK.e] = -^{/i,k«} = °! 

where the bracket {•, •} denotes the Poisson bracket associated to the symplectic form u. Due to 
this equality, we can write (see for instance [ AMR99| , Corollary 4.1.27]) the following expression 



for Gt, the flow of ^_k« : 

G t (v) = lim (F t/n o K_ t/n ) n {v) = (K^ t o F t )(v) = exp • F t (t;). 

n— >oo ' ' 

Since by hypothesis the point v is periodic for Gt with period r, we have that 

v = exp — t£ • F T (v), 

or, equivalently, 

F T (v) = exp r£ • v, 
12 



as required. ■ 

Using the previous proposition, we will reduce the search for RPOs of a generic one- 
parameter family of G-Hamiltonian systems (V,u, h\) that satisfies conditions (HI), (H2), (H3), 
and (H4), to the search for periodic orbits of the vector fields of the form X hx _^, and will 
prove the following result: 

Theorem 3.3 Let (V,u>,h\) be a one-parameter family of G -Hamiltonian systems that satis- 
fies conditions (HI), (H2), (H3), and (H4). Suppose that G contains a Lie subgroup H of 
positive dimension with associated equivariant momentum map K : V — * h*. Let U Uo be the res- 
onance space with primitive period T Uo . Then, for each £ € f) whose norm ||£|| is small enough, 
there are at least, in each energy level nearby zero and for each value of the parameter A near 
A , as many relative periodic orbits as the number of equilibria of a G^ x S 1 -equivariant vector 
field defined on the unit sphere on Vq. The symbol G^ denotes the adjoint isotropy subgroup of 
the element £ € h, that is, 

G^ = {geG\ Ad 9 £ = £}. 



Remark 3.4 If we are just interested in looking for purely periodic orbits it suffices to use 
Theorem 3J3 with £ = 0. Conversely, if we use this result with a value of the parameter £ ^ 
we cannot conclude that the predicted RPOs are not trivial, that is, that they are not just 
periodic orbits. This point will become much clearer in the examples presented in the following 
sections. ♦ 



Remark 3.5 In terms of practical applications, the relevance of Theorem |3.3| is given by the 
fact that the estimate that it provides in terms of the number of equilibria of an equivariant 
vector field on the sphere can sometimes be calculated via topological arguments, as we will 
see later on. ♦ 

Proof We will work in the basis of the resonance space U Vo provided by the equivariant 
Williamson normal form, in particular we will use the matricial expressions fl2.3| ), which are 
consistent with the decomposition U Uo = Vq®V\ presented in ( 2.1l|) . Recall that the subspaces 



Vq and V\ are G x S^-invariant. Abusing the notation a little bit we will use the symbol £ to 
denote both an element of the Lie algebra f) C and its representation on Vo and V\. Using 
Lemma [O] we can write, for each v = vq + v i € U Uo represented in the previously mentioned 
basis, 



t 0\(v 
{JU 



Note that, also by Lemma EO, that the matrix £ is skew-symmetric, £ T = — £, therefore normal 



and hence diagonalizable. The same Lemma implies that the linear map £ : Vq — > Vq associated 
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to £ € f) commutes with ^2n, [£> J2n] = 0, and consequently these two endomorphisms can be 
simultaneously diagonalized. 
We recall that, 

(K(v),0 = \u{t-v,v) = l(v , Vl ) ( ° J ) ( £ 

In particular, 

d 2 K«(0) = 

We start the proof by defining the R x fj— parameter family of Hamiltonian functions given 



£ 

< o 



by 

hxt = h x - K«. 

Due to the hypotheses on the family h\, the quadratic nature of the momentum map K, 
and the fact that h\o = h\, the family h\ £ satisfies the hypotheses of the Normal Form 
Reduction Theorem [ VvdM9q] . Therefore, a new family can be constructed such that, for 
any value (A, £) of the parameters, the Hamiltonian h\£ is if? 1 — invariant with respect to the 
action generated by the semisimple part of the linearization at zero of X/ lXo = , that is, 

(9,v) i ^ e^ A °v, 9 e S l , with 

AS _ f "ohn \ , . 

The Normal Form Reduction Theorem guarantees that the 5 1 -relative equilibria of h\£ are in 
correspondence with the periodic orbits h\£ which, by Proposition 3.2, are RPOs of h\. The 
quadratic nature of the momentum map K and its S 1 -invar iance imply that h\£ can be chosen 
to be of the form 

h\£ = h x - K 5 , 

with h\ the normal form for the family h\. 

As a result of thesejaremises, the RPOs that we are looking for will be given by the critical 
points of the function h\ — — 3^ +a , that is, the elements (v, a, A, £) € U Va x R x R x t) for 
which the function 

FHv, a, A) := V^ o (hx - - J< +a ) (v) (3.2) 



has a zero. As customary, the gradient in the previous expression is constructed using the inner 
product introduced in Lemma [2.1|. 
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Lyapunov Schmidt reduction and the bifurcation equation The linearization L*> : 
U Va —> U Va of the equation (^) at the point (0,0, A o ,0) produces, in the usual basis, the 
expression: 

LC = d 2 (^-jC)(0) = d 2 (/ iA -^)(0)=f ° (1 "; KJ 4 (3.3) 

By looking at this matricial expression we see that it is possible to Lyapunov-Schmidt reduce 
the bifurcation problem posed in (3.2) whenever £ = v Q , which we will assume in the sequel. 
In those circumstances ker L u ° = Vq, Iml"° = V\. Let P : U Uo —> Vq be the G x S^-equivariant 
projection associated to the splitting U Vo = Vb©Vi. The equation (I — ¥)F U ° (vq + v i, a, A, £) = 

(I — P)V[/„ o ^/ia — — J Uo+a ^j (vq + vi) =0 defines, via the Implicit Function Theorem, a 
function v\ : Vq x R x R x f) — > Pi, such that 

(I-P)i^ (uo + vi(«o,a,A,e),a,A,e) = $-V)Vu Vo (£ - & - J^+j K + ui(«o,a, A,£)) = 0. 

(3.4) 

Notice that the function v\(vq, a, A, £) is x S^-equivariant, since this is the symmetry under 
which F v ° is equivariant, that is, for any g E G^ x S 1 , we have that v\(g ■ Vq,ol, A,£) = 
g • ui(u ,a, A,0- 

The final Lyapunov-Schmidt x S^-equivariant reduced bifurcation equation, whose 
zeros provide us with the RPOs that we are after, is given by-B:VoxRxRxf) — > Vq, where 

B{v ,a,X,0 = FF»°(vo + v 1 (v ,a,\,0,a,\,0=WV u „ a (h~\ - K« - r° +a ) (vq + vi(v ,a,X, 

= Vu^lhx-K^-r^ivQ + v^a^,^)) (by (3.4)). (3 

We collect the main properties of the reduced bifurcation equation in the following 

Lemma 3.6 The reduced bifurcation equation ( |3.5| ) is G^ x S 1 -equivariant with respect to the 
action of this Lie group on Vq and it is the gradient of a G^ x S 1 -invariant function defined 
on Vq, that is, 

B(v , a, A, f) = V Vo g(v , a, A, £), 

where the function g : Vq x Lie(S' 1 ) x R x f) — > Vo is defined by 

g(v ,a,X,0 = (hx- .F° +a -K%o + ^o,a,A,£)). 

Proof The G^ x 5 1 -equivariance is a direct consequence of the construction of B. As to the 
gradient character of B, note first that for any w E V\ we have that 

{F v ° {v + v 1 (v ,a,\,S),a,\,Z),w) = (F v ° (v + v x (v , a, A, £) , a, A, f ) , (I - P)u>) 

= ((I-P)F'"'(«o + u 1 (« ,a,A,e),a s A,0,«>> =0 (3.6) 
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where the last equality follows from the construction of the function v\ through expression (| 
Now, let u £ Vq arbitrary. We write: 

(B(v ,a,X,O,u) = (FF u °(v + vi(v ,a,\,0,a,\,0,u) 
= (F Uo (v + v 1 (v ,a,\,£),a,\,Q,u) 

= (F Uo (v + v 1 (v ,a,\,£),a,\,£),u + Dv^ivo,!*^^) ■ u) (by (^)) 
= <Vtf„ (hx - J 1/0+0 - K*)(u + vtivo, a, A, £))>« + Ovb«i(«o» <*> A - ■ u) 
= d(hx - J 1/0+0 - K*)(u + wi(«o, a, A, 0) • (« + D VoVl {vo, a, A, ■ u) 
= dg(v ,a,\,£) ■ u = (V Vo g(v , a, A, £), u), 



as required. This construction is a particular case of the one carried out in [GMSD95] and [CLOR99] . 
T 

Notational simplification: In order to make notation a little bit lighter we will assume in 
the rest of the proof, without loss of generality, that the system has been scaled in such a way 
that u = 1 and A = 0. 

The following lemmas provide a local description of the reduced bifurcation equation that will 
be much needed. 



Lemma 3.7 The function v\ introduced in ( |3.4j) has the following two properties: 

(i) vi(0,a,\,£) = for all a, A G M, and £ G fj. (3.7) 

(ii) £V ui(0,a,A,0 = ~^hn ~ (1 + a) ~^ (A) Jan - (3-8) 

p(A) p(A) p(A) 

Proof Part (i) is a consequence of the uniqueness of the solutions provided by the Implicit 
Function Theorem. The proof of part (ii) is supplied in the Appendix, Section |5.2| . T 

The proof of the following lemma is a lengthy but straightforward computation. 

Lemma 3.8 Let B(vo,a,A,£) be the reduced bifurcation equation, then: 

(i) 

n B , n x ^ a(X)p(X)-r i (X)-((l + a)-iP(X))\ 2 [(1 + a) - ip(X)] £ 2 
D Vo B(0, a, X, £J = -— hn H ttt hnt, + 



p{X) ™ p(A) ™ p(A)" 

(ii) T/ie principal part of the reduced bifurcation equation is given by the expression: 

B(vo, a, A, f ) = (A(t'(0) + a> - f 2 «o " 2aJf 2 n£vo - 2^ / (0)aAw 

+ 2/(0)AJ 2 n^o + C (4 3) ) + ^-o.t, (3.9) 

where C (v^j is the trilinear operator obtained by taking the gradient of the fourth order 



term in the vo~expansion of h\ a (vq + t>i(t>o, 0, 0, 0)). 
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We now write the reduced bifurcation equation in polar coordinates, that is, we define 

B p (r,Uo,a,\,£) = B(ru ,a, A,0, 
where r£l and uq £ gd 1111 ^!}-^ "We introduce the function 

F(r,uo,a,\,£) = . 

r 

By looking at ( |3.9| ) it is clear that the function F is smooth at the origin, F(0, 0, 0, 0, 0) =0 
and that D x F(0, 0, 0, 0, 0) = cr'(0) / 0, by hypothesis (H4). Therefore, the Implicit Function 
Theorem guarantees the existence of a smooth function A(r, uo, a, such that A(0, 0, 0, 0) = 
and F(r, uo,a, A(r, uo,a, 0,0 = 0. This equality implies that if we substitute the function 
A(r, Uq,oi,£) on the reduced bifurcation equation, this time considered as a vector field on Vb, 
we obtain a new (a, ^-parameter dependent vector field 

G(r, u , a, = B 6 (r, u , a, A(r, u , a, 0,0 (3-10) 

which due to the fact that {Bf,(r, uq, a, A(r, no, a, 0, 0, u o) = is, for each small enough fixed 
value of r, a x S^-equivariant vector field on the sphere on Vb of radius r, whose zeroes 
constitute solutions of the reduced bifurcation equation. ■ 



Method for the optimal use of Theorem 3.3 The optimal and most organized way to 
apply Theorem 3J3 consists of using the estimate it provides in the fixed point subspaces Vq 



corresponding to the various subgroups H in the lattice of isotropy subgroups of the G x S 1 - 
action on Vb, replacing the group G x S 1 by N(H), which is a group that acts on (not 
necessarily in an irreducible manner). The symbol N(H) denotes the normalizer of H in G x S 1 
and Vff 1 is the vector subspace of Vo formed by the vectors fixed by H. We make more explicit 
this comment in the following paragraphs. 

Let H be a subgroup of G x S 1 . If n : G x S 1 — > G denotes the canonical projection and 
tt(H) =: K C G, Proposition 7.2 in [GSS88| guarantees the existence of a group homomorphism 



6 : K -> S 1 such that 

H = {(k, 9(k)) eG x S 1 \k £ K}. (3.11) 

In our discussion we will be concerned with spatiotemporal symmetries, that is, subgroups 
H of G x S 1 for which the homomorphism 9 : K — > S 1 is nontrivial. Using the characteriza- 
tion ( |3.11| ) it is straightforward to see that 

AT (if) = N G (K) x S 1 . 

The Ng {K )-action on U„ is globally Hamiltonian with momentum map : Lie (Ng(K)Y 

given by the restriction of the G-momentum map to U^, that is, for any v G and any 
^ E Lie (Nq(K)), we have that 

(K»,0 = (K(n),0- 

The same statement applies to the S^-action. Using these objects we can reformulate Theo- 



rem 3^ on the fixed point spaces Vq. 
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Corollary 3.9 Let (V,u,h\) be a one-parameter family of G-Hamiltonian systems that sat- 
isfies conditions (HI), (H2), (H3), and (H4). Let H be a spatiotemporal isotropy sub- 
group of the G x S 1 -action on Vq, such that dimV^ = 2k and K := n(H). Then, for each 
£ G Lie (-/Vg(-?0) whose norm ||£|| is small enough, there are at least in each energy level nearby 
zero and for each value of the parameter A near X Q , as many relative periodic orbits as the 
number of equilibria of a Ng(K)^ x S 1 -equivariant vector field on the unit sphere on V H . The 
relative periods of these RPOs are close to T Uo , and their phase shifts are close to expT Uo ^. 
The symbol Nq(K)^ denotes the adjoint isotropy subgroup of the element £ £ Lie (Ng(K)), 
that is, 

N G (K)t = {ge N G (K) | Ad 3 e = £}• 
The mapping tt : G x S 1 — > G denotes the canonical projection. 



As we already said, both the previous result and Theorem 3.3 can be used to look for purely 
periodic motions by taking in their respective statements £ = 0. There is a situation of special 
interest: 

Periodic orbits with maximal isotropy subgroup: Let H be a maximal isotropy subgroup 
of the G x S^-action on Vq. In the presence of maximality we have at our disposal the following 
convenient result: 

Lemma 3.10 Let H be a maximal isotropy subgroup of the compact G x S 1 -action on Vq. Let 
N be the Lie group N(H)/H and N° be the connected component of the identity of N. Then 
either 



(i) ~ S l , and N/N° = {Id} or N/N° ~ Z 2 , or 

(ii) N° ~ SU{2) and N ~ SU(2). 

In the first case we say that H is a maximal complex subgroup. In the second case we say that 
H a maximal quaternionic subgroup. 

Proof It is a straightforward combination of the general result for linear actions of compact Lie 



groups [ Bre72 , G83 , GSS88 | with Proposition 12.5 in [ GoSt85 | that eliminates the possibility 
of having real maximal isotropy subgroups when the compact group in question is G x S 1 . ■ 

Using the previous lemma and an additional genericity hypothesis, the estimate given in 



Theorem 3.3 can be made very explicit: 



Corollary 3.11 Let (V, u, h\) be a generic one-parameter family of G-Hamiltonian systems 
that satisfies conditions (HI), (H2), (H3), and (H4). Let H be a maximal isotropy subgroup 
of the G x S 1 -action on Vq such that dim(V^) = I 0. Then: 
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(i) If N° ~ S 1 there are at least 1/2 (if N/N° = {Id}) or Z/4 (if N/N° ~ Z 2 ) branches of 

periodic solutions with isotropy H coming out of the origin as one varies the parameter 
X, with periods close to T Uo . 

(ii) If N° ~ SU(2) there are at least I /4 branches of periodic solutions with isotropy H coming 

out of the origin as one varies the parameter X, with periods close to T Uo . 



Proof We will adapt to our problem the approach followed in [ CKM95| , Koe95| ] for rotating 



waves. The main idea behind the proof consists of using the maximality hypothesis to give 
a numerical evaluation of the estimate in Corollary |3.9| , that is, the number of equilibria of a 
7V(.£f)/ii-equivariant vector field on the sphere 

More especifically, let G H := G\ v b be the restriction of the vector field G on Vq, defined 



in ( 3.10 ), to the fixed point set V H , and G^ 1 (uq, a) := G H (r, uq, a) be the iV-equivariant vector 
field on S 1 ^ 1 obtained by fixing r in the mapping G H (note that in our case £ = since we are 
looking for periodic orbits). The zeroes of this vector field are in one to one correspondence 
with the solutions that we search. Due to the maximality hypothesis on the subgroup H, the 
iV-action on the sphere is free and therefore the corresponding orbit space S l ~ l /N is a 
smooth manifold onto which we can project the iV-equivariant vector field G^. Let G^ be the 
projected vector field. Due to the genericity hypothesis in the statement, the Poincare-Hopf 
Theorem allows us to say that G^ has at least xiS 1 ^ 1 /N) equilibria, where x denotes the 
Euler characteristic. These zeroes lift to equilibria of the restriction of the reduced bifurcation 



equation to V H , due to the gradient character (see Lemma 3J:) of B and consequently of its 



restriction to V H . 

In order to conclude our argument it is enough to show that xiS 1 ' 1 / N) corresponds to 
the estimates provided in the statement of the theorem. In the first case, when iV ~ S , the 
dimension of V H is necessarily even (we will write I = 2k for certain k € N) and there are two 
possibilities: the quotient N/N° is either {Id} or it is isomorphic to Z 2 . If N/N° = {Id}: 

X (S l - l /N) = x((S l - l /N°)/(N/N )) = xiS^/S 1 ) = x^ 1 ) = k = l -. 

If N/N° ~ Z 2 : 

xiS^/N) = xdS^/N^/iN/N )) = X (CP fc -Vz 2 ) = \ = \, 

where we used the well-known fact that if G is a finite group acting freely on a manifold M, 
then (see for instance | Kaw9l| , Corollary 5.22]) 



X 



(M\ _ X (M) 



\G J \G\ ' 

Finally, if H is maximal quaternionic then / = dim(V H ) = Ak for some k £ N, necessarily, and 
xiS^/N) = X (S 4k - 1 /SU(2)) = xOHIP^ 1 ) = k = ~ 
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The calculation of the Euler characteristic x(HP fc l ) of the quaternionic projective space is 
made using an argument based the spectral series of Leray (see for instance | BT82f ). 



The computations that we just carried out give us periodic orbits for a fixed r. Moving 
smoothly this parameter we obtain the branches required in the statement of the theorem. ■ 

4 Bifurcation of non-periodic relative periodic orbits in the 
presence of extra hypotheses 



The tools presented in Theorem |3.3| for the search of RPOs based on topological methods 
produce estimates that, as we will see in the following examples, have some limitations, in 
particular we have no example where it guarantees the bifurcation of non periodic RPOs. This 
circumstance has motivated us to use a more analytical approach under dimensional hypotheses 
that are satisfied in very relevant situations. A detailed study of the bifurcation equation in the 
presence of these hypotheses will provide us with sharper estimates that completely describe 
all the bifurcation phenomena that we see in the examples. 

4.1 Motivating example: two coupled harmonic oscillators subjected to a 
magnetic field as an example of symmetric Hamiltonian Hopf bifurcation 

We consider the system formed by two identical particles with unit charge in the plane, sub- 
jected to identical harmonic forces, to a homogeneous magnetic field perpendicular in direction 
to the plane of motion, and to an interaction potential that will preserve certain group of sym- 
metry. We will denote by (q%, (72) the coordinates of the configuration space of the first particle 
and by (53, those of the second one. If 7 is a constant that determines the intensity of the 
magnetic field, it is easy to see that the Hamiltonian function of the system described above is 

H (q ' p) = + ^ + p z + ^ + ( -hi ~ I) ^ + q * + q * + ^ 



7 7 ... 

H (pm - Pzqi) H (f>3<?4 -P4Q3) + fiA^l^D, (4.1) 

m m 

where 

A = Qi + Qi+2, A = Pi + Pi+2, ^3 = PM+2 ~ Pi+2Qi, A = qiPi + qi+2Pi+2, i G {1, 2}, 

and / is a higher order function on its variables that expresses a non linear interaction between 
the two particles. 

This system has, for all values of the parameters 7 and k, an equilibrium at the point 
(<Zi> <&> <72) Q&iPiiP3iP2,P4) = (0,0). The linearization of the dynamics at that point is repre- 
sented by the matrix (the coordinates are ordered as in the previous equality) 



_7J J_j \ 

1. -a* ■ (4 ' 2) 
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whose eigenvalues are 



1 



Afc = ± — \J km — 27 2 ± 27 \/ 7 2 — km. 

If we move the parameter k around the value k a = 7 2 /m these eigenvalues present a Hamilto- 
nian Hopf behavior like the one depicted in Figure 1.1. 

We now study the symmetries of the system. Note that after the assumptions on the 
interaction function /, the system is invariant under the canonical S^-action given by the 
lifted action to the phase space of 



/ 

V 



cos (p — sin ip 
sin (p cos <p 









\ 



COS If 

sin ip 



— sin ip 
cos tp 



where q = (qi,q3,q2, qi), and by the transformation 



f qi\ 

q-i 
93 

V 94 / 



/ ft \ 

92 
-93 
V "94 J 



The momentum map K : M 8 — > R associated to the 5' 1 -action is given by the expression 

K(q, p) = £> 3 9i - 93Pi - £>2<74 + £>4<72- 

If we now look at the linearization (|4.2j ) evaluated at the Hopf value k = 7 2 /m we see 
that in this case Vq consists of the points of the form (qi, q%, q 2 , <?4, 0). The 5 1 -action on Vq 
generated by the semisimple part of Ak a can be written as 



7J4 



In order to better study the group actions on Vo we will perform a linear change of variables. 
Let (z±, Z2, z~i, Z2) be the new (complex) coordinates, given by 



zi = qi + 94 + iq2 - iqz 

Z2 = qi - qi + iq2 + iqz 

zi = qi+qi- iq2 + % 

z~2 = qi - qi - iq2 - iqz- 



(4.3) 



If we take as new angles ipi and ^2, defined by: 

7 n 



1p2 = <f 



m m 

we realize that the previously introduced actions form a 0(2) x S^-action on Vo that takes the 
following convenient simple expression: 



(Vi,V 2 )-(^ 2 ) = (e #1 



zi,e^z 2 ) 



and 



r ■ (zi,z 2 ) = (z2,zi). 



(4.4) 
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That is, we have shown that the system of two coupled harmonic oscillators subjected to a 



magnetic field, whose Hamiltonian is given by (iA) can be taken as an example of Hamiltonian 
Hopf bifurcation with 0(2) x S 11 — symmetry, which we will study in full generality in the following 
subsection. 



4.2 RPOs in Hamiltonian Hopf bifurcation with 0(2)— symmetry 

Having the example in the previous section as a motivation we will study in what follows the 
RPOs that appear in a Hamiltonian Hopf bifurcation phenomenon in the presence of a 0(2)- 
symmetry. The simplicity of this symmetry will allow us to explicitly write down the principal 
part of the reduced bifurcation equation in full generality, and to read off directly from it the 
RPOs that we are looking for. 

We start by recalling that in the canonical coordinates introduced in (|2.3[) the principal 



part of the reduced bifurcation equation is, by Lemma 3.8, equal in our case to 

B(vq, a, A, £) = (A<t'(A ) + a 2 z^)u - t, 2 v - 2av J^v G - 2^' (\ o )v o a\v 

+ 2V / (A )AJ 4 ^o + IPd 4 /i Ao (0) (v^ ] ) i!u,l.. 



We now rewrite this expression in the coordinates (z±, z 2 , z\, z 2 ) in which the 0(2) x S^-action 
looks like ( f4.4[ ). In doing so we need to express in these new coordinates the matrix J^, which 
can be easily achieved by using the explicit expression of the change of variables (|4.3| ). Indeed, 
we have that in those coordinates 



J 4 £ = 



/ 1 

o 



V o 



\ 



•i / 



and therefore, the first terms of the expansion of the reduced bifurcation equation are: 

-Z2 ) 



B(z,a,X,0 = (Aa'(A ) + aX) 



Z 2 
V *2 J 



2v Q ai 



\ 



+ 



( (a\ Zl \ 2 + b\z 2 \ 2 )z 1 ^ 
(a\z 2 \ 2 + b\z 1 \ 2 )z 2 
(a|zi| 2 + b\z 2 \ 2 )zi 
V (a\z 2 \ 2 + b\zi\ 2 )z- 2 ) 



+ ... 



(4.5) 



where the coefficients a and b are related to the fourth order terms in the expansion of the 
Hamiltonian, that is, Pd 4 /i Ao (0) (^o^) ■ ^ n or der to keep the simplicity of the exposition we 
will assume that these two coefficients are non zero and non equal (otherwise we would have 
to go to higher orders in expression ( f4.5p ). The RPOs that we are looking for are given by the 
solutions of the system of equations: 

= (Xa'(X ) + n 2 v 2 )zi - 2v Q aiz x + (a\zi\ 2 + b\z 2 \ 2 )zi + ... (4.6) 
= (Ac/(Ao) + a 2 u 2 )z 2 + 2v Q aiz 2 + (a\z 2 \ 2 + b\ Zl \ 2 )z 2 + ... (4.7) 
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Since <t'(A ) 7^ 0, equation ( |4.6| ) can be easily solved by dividing the expression by Z\ and then 
using the Implicit Function Theorem to define a function 



cr'(A 



[-a 2 v 2 + 2i/ a^ - (a|zi| 2 + b\z 2 \ 2 ) + ■■■'_ 



(4.8) 



that substituted into ( |4.6| ) solves it. Hence, plugging (4.8) into ( |4.7|) we reduce the problem 
to solving a scalar equation which can be done again via the Implicit Function Theorem: we 
divide ( |4.6|) by Z\ and Q4.7|) by Z2 and obtain: 



= 4z; a£ + (a - 6)(|z 2 | 2 - ki| 2 ) + ■ • • 



(4.9) 



Since for equivariance reasons z\ and z 2 always appear in the tail of the previous expression 
as combinations of |zi| 2 and | ^2 ! 2 ; the hypotheses on the coefficients a and b allow us to solve 
this final scalar equation by defining a function 



\z 2 \ 2 = W? (\zi\ 2 ,a,(,) = \zx\ 



4f ai£ 



(4.10) 



As we illustrate in Figure 4.1, the solution (4.10|) predicts, for each fixed value of the norm 



|z 2 | 2 






\^" lz,l! 




/ '|z,| 2 




/\ \z,\ 2 -(4vj(a-b))af 



Figure 4.1: Parameterization of the RPOs in the Hamiltonian Hopf bifurcation with 0(2)- 
symmetry. 



M 2 + 



z 2 \ 2 a one-parameter family of RPOs that are obtained by making vary the product 



a£. More explicitly, and using Figure 4.1, suppose that 



kil 2 + 



z 2 



is fixed, then, the intersection of the lines 
with the abovementioned RPO. 



z 2 



1 — \z\\ and \z 2 \ 



1 and that the value a£ 
\Z\\ — — provides us 
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Notice that all these RPOs cannot be predicted by merely using Theorem 3.3 since even 
though all the hypotheses needed in this result are fulfilled, it only predicts two RPOs for each 
value of the norm \zi\ 2 + \z2\ 2 - Moreover, one cannot be sure that these are not just periodic 
motions since the predicted orbits could lie in fixed spaces of maximal isotropy, thereby implying 
their periodicity. 

Remark 4.1 In contrast with the Hamiltonian case, the Hopf bifurcation of nontrivial RPOs in 
the dissipative case with 0(2) symmetry is subjected to the presence of a cubic order degeneracy 
in the normal form. Unfolding this singularity leads to a codimension two bifurcation problem 
where the RPOs appear as a secondary branching from the primary branches of periodic orbits. 
♦ 



4.3 Hamiltonian Hopf bifurcation of RPOs for reduced integrable systems 

The analysis performed in the previous section dealing with the Hamiltonian Hopf bifurcation 
with 0(2)-symmetry is a particular case of a more general situation. Indeed, the main feature in 
that example that allowed us to carry out a by hand in-depth study of the reduced bifurcation 
equation was the coincidence of one half the dimension of the reduced space Vo with the 
dimension of the symmetry group 0(2) x S . We will see in this section that whenever we 
are in the presence of this reduced integrability hypothesis an analysis in the same style can be 
performed. 

More explicitly, all along this section we will be dealing with (V,u,h\), a one-parameter 
family of G-Hamiltonian systems that satisfies conditions (HI), (H2), (H3), and (H4) such 
that if 4n is the dimension of the resonance space U Vo with primitive period T Va , then the rank 
of G x S 1 equals n, that is, the maximal tori of the Lie group G x S 1 have all dimension equal 
to n. 

Let T n_1 C G be a maximal torus of G, and let £ £ t™" 1 be an element in the Lie algebra 
of Tpn-i. Like in the previous section we can find coordinates in which the action of T n 1 x S 1 
looks simple. Namely, there exists a set of complex coordinates (z\, . . . ,z n ,z\,... , z n ) (and 
conjugates) for Vq and a set of angular coordinates (£i, . . . , £ n -i) for the torus T n_1 , for which 
the T n_1 -action looks like 

(e* 1 ,... ,e^- 1 )-(z 1 ,... ,z n ) = (e* 1 z 1 ,... , e^-^ n _i, e < c ^+~ +c "-^-^ z n ) , 

where the coefficients c\, . . . , c n _i are rational constants. If we incorporate the S^-action using 
these complex coordinates, the T n_1 x S^-action looks like 

{e ixi \... ,e^-\e iQ ) • (zi,... ,z n ) = (e i($1+Q) zi, . . . , e i( *"- 1+a) * n _i, Zn ). 
Let us now set 

1>j = & + » j = l,-..,n-l (4.11) 

V>n = Cifi + • • • + Cn-iU-l + OL. (4.12) 
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Under the condition 



ci H h c n _i / 1 



(4.13) 



these relations define a change of coordinates on the n-dimensional torus T n 1 X S , and in 
these new coordinates the action can now be written in the very simple fashion 



(e^ 1 ,... ,e**») -fa,... ,z n ) = (e^ Zl ,... ,e^z n ). 



J4>n 



(4.14) 



Notice that under Condition ( |4. 13| ) the ring of invariant polynomials for this action on Vq is 
generated by the quadratic invariants ttj = ZjZj, j = 1, . . . ,n, and that the strata of this action 
are obtained by setting some of the Zj's equal to while keeping the others different from 0. The 
orbit space for this action can be identified with the positive cone in R n {(tti, . . . ,ir n ) j ttj > 
0, j = 1, . . . ,n}. 

Recall now that in the canonical coordinates introduced in ( |2.3[) , the principal part of the 
reduced bifurcation equation is, by Lemma 3.8, equal to 

B(v , a, A, £) = (A(j'(A ) + a 2 iy 2 )v - £, 2 v - 2a^ o J 2 n^0 - 2ip' (\ )u aXvo 

+ 2^'(Ao)AJ 2 ^o + C UfA + h.o.t., (4.15) 

From ( 4.14| ) it is clear that the matrices J 2 n ; $2n£, and £ 2 in ( 4.15| ) take, in these newly 
introduced coordinates, the form: 





f 1 






§2n — 






hnC = 











( -fa 



V 



n-l 



\ 



-i>n J 



\ 



o \ 



^ J 



Using these new coordinates and the T n 1 x S 1 equivariance of B, we rewrite the components 
of ( [4.15| ) (we omit the complex conjugate part). For i G {1, . . . , n} we have: 



Bi(z,a,\,0 = Aa'(A ) + $ + Q (M 2 ,... , \z n \ 2 ) + h.o.t. z 



where 



Cil | ^1 1 ~i~ ■ ■ ■ Cin I Zn I ) 



Cij £ 



i,j G {1, . . . ,n}, 



We can now state a theorem about the bifurcation of RPOs. 
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Theorem 4.2 Let (V,uj,h\) be a one-parameter family of G-Hamiltonian systems that satis- 
fies conditions (Hl) ; (H2), (H3), and (H4). Suppose that: (i) the dimension of Vq equals 
twice the rank n of G x S 1 ; (ii) the condition ( \j-lS ) on the torus action is satisfied. Then, if 
the matrix 

A = (c^j — Cij) , 1 < i < n, 1 < j < n — 1, 

has maximal rank n — 1, there exists a family of RPOs with n different frequencies which 
bifurcates from the trivial solution as X crosses Xq. 

Remark 4.3 The condition on the matrix A is not generic, because the values of the coeffi- 
cients c^j are constrained by the G-equivariance of the operator C. 

Proof Since we are looking for solutions with Z{ ^ 0, for all i, we can factor out Z{ in each 
equation Bi = 0. The resulting equations read, for i £ {1, . . . , n}, 

= \a'(X ) + i;f +Cl(\ Zl \ 2 ,... ,\z n \ 2 ) + h.o.t. 



These equations are simply those that we would have obtained by projecting first ( 4.15 ) on 
the orbit space corresponding to the toral action; in all that follows we will denote \zi\ 2 by 7Tj. 
Since by hypothesis (H4), a'(X ) ^ 0, we can solve any one of these equations for A. Let us do 
so for the equation with i = n. By substituting A by the resulting expression in the remaining 
equations, we have reduced the problem to solving a system of n — 1 equations which, at leading 
order, have the form 

= ip 2 - ipl + (t5i - Q)vri H h {cTn - c^ n )ir n + h.o.t.. (4.16) 

If the matrix A has maximal rank, we obtain a unique family of solutions of the system (|4.16[) , 
for which n — 1 of the 7Tj's depend smoothly on the remaining one and on the parameters tpj, 
j = 1, . . . ,n. In order to fix notations, let us assume without loss of generality that we have 
obtained 7Tj = TTifyi, . . . ,ip n ,n n ) for i = 1, . . . ,n — 1. These solutions still have to lie inside 
the orbit space, that is, we still have to check the additional conditions 7Tj > 0. However, since 
the ipj's are free parameters, the quantities ip 2 — ip 2 can take any real value. Therefore, if 
we set z n = 0, we can always find values for the ^-'s such that 7Tj > for i = 1, . .. ,n — 1. 
These inequalities are still satisfied if the f/'j's are close enough to these values and 7r n > is 
close enough to 0. Finally, since these solutions lie on the principal stratum for the action of 
T n_1 x S , the corresponding RPOs have n different frequencies which depend smoothly on A. 
■ 

Remark 4.4 In the problem with 0{2) symmetry analyzed in Section 4.2, we have n = 2, 



c\ = —1, en = C22 and c\i = C2i- The hypotheses of Theorem [L2| are therefore generically 
satisfied in this case. ♦ 



Remark 4.5 As it was already the case with Theorem |3.3| , Theorem |4.2| still applies if instead 
of G x S 1 acting in Vq, we consider the group N{H)/H acting in V H , where H is an isotropy 
subgroup of the G x S^-action. In the next section we shall see an application of this remark. 
♦ 
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4.4 Hamiltonian Hopf bifurcation with SO (3) symmetry 

Hopf bifurcation problems with SO(3) symmetry for dissipative systems have been investigated 
by several authors in the case in which the eigenspaces associated with the critical eigenvalues 
is the direct sum of twice the five dimensional (real) irreducible representation of SO (3) (see 
| GoSt8EH ), |oRo89H , [|MRS88|1 , and |Le97[| ). This is the simplest possible case with SO (3) 
symmetry which does not reduce to Hopf bifurcation with either trivial or 0(2) symmetry. 
Nevertheless, it leads to a normal form in a ten dimensional real vector space. The list of 
solutions with maximal isotropy, hence purely periodic ones, has been given in |GoSt85 ] and 
in [MRS88]. However, the most interesting feature of this problem is the possibility of having 
a bifurcated branch of RPOs in a six dimensional subspace. This was first found by [|IoRo89|l . 
Another approach was taken by [ Le97f| (using orbit space reduction) who did not recover the 



result of [IoRo89|. This remark shows the level of difficulty found in obtaining direct branching 
of RPOs via Hopf bifurcation for equivariant vector fields. In the Hamiltonian context, a 
related work by Haaf, Roberts and Stewart [HRS92] has shown the existence of families of 
periodic orbits with maximal isotropy for a Hamiltonian in M 10 which is invariant under the 
same SO(3)-action. 

In the sequel we investigate the Hamiltonian Hopf bifurcation with 50(3) symmetry, when 
the subspaces Vq and V\ are associated with this ten dimensional representation and we shall 
see that, in this case, Theorem [1^ applies and shows the existence of several families of RPOs. 

Let 

V = V © Vx ^ M 10 . 

We identify M 10 with M 5 <8> C and consider the action of SO (3) on R 5 given by its irreducible 
representation on the space of spherical harmonics of degree 2. Equivalently, we may identify 
M 5 with the space W of 3 x 3 real symmetric matrices with trace 0, and consider the action of 
SO (3) on W defined by 



p A (M) = A' 1 MA, A e SO(3), M €W. 



This definition extends naturally to 



C with the same formula, M now having complex 



coefficients. We shall therefore identify in all that follows Vq with W <8> C. The S -action on 



Vo is simply defined as multiplication by e in C, that is, ■ M := e M. 
Any M G W <g> C decomposes uniquely as 



M 
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where 



Bo 



1 

10 

2 

1 i 
i -1 

























I— 1 


i 


o / 



Bi 



B. 



Bo. 



(4.17) 
(4.18) 



We now list the isotropy types for the action of G = SO (3) x 5 1 on Vq that we have just 
defined. We use the presentation and results of | HRS92(| . Figure [T^ shows the isotropy lattice 
of the G-action and the dimension of the corresponding fixed-point subspaces. 

dim V H 




0(2) D 4 T SO(2) <1) SO(2) (2) 2 




6 
10 



Figure 4.2: Isotropy lattice of the 50(3) x 5 1 -action. 



Notations: if is a subgroup isomorphic to H C 0(3) but such that H H 5 1 7^ 1 (here 1 is the 
trivial group). In particular, Z ra is the group generated by (R n ,— 27r/n) G 50(3) x 5 , where 
R n is a rotation of angle 2ir/n. 

By Corollary [T^, for any subgroup H with dimiVQ 1 ) = 2 there exists a branch of periodic 
solutions having this symmetry. Let us now consider those isotropy subgroups having a fixed- 
point subspace of dimension 4. The largest subgroup acting faithfully in Vq is N(H)/H. We 
list below the different N(H)/H for the non maximal isotropy subgroups: 

1. N(D 2 )/D 2 ~D 3 x 5 1 ; 

2. N(Z 4 )/Z 4 ~ N(Z 2 )/Z 2 ~ 0(2) x 5 1 ; 

3. iV(Z 3 )/Z 3 ~ 50(2) x 5 1 ; 

4. iV(Z 2 )/Z 2 ~ 0(2) x 5 1 ; 
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5. JV(1) = G. 



In Case 1 we see that solutions with isotropy D 2 are always periodic. Case 2 corresponds to 
the problem described in Section 4.2 (Hopf bifurcation with 0(2) symmetry). It was noticed 
in fHRS92| 1 that the equations in V Q H do not degenerate despite the fact that they come from 
a system with higher symmetry, which leads us to conclude that families of RPOs with two 
frequencies and with spatio-temporal symmetry Z4 as well as r L 2 do generically bifurcate. Case 
3 falls in the framework of Section 4.3, since the symmetry group is 50(2) x S . However, 

because the equations in Vq 3 are the restriction in that subspace of a system with higher 
symmetry in Vq, we need to compute the cubic order terms in order to insure that no "hidden" 
degeneracy occurs. We use the argument proved in [HRS92]; we can choose Z3 so that, by 
introducing complex coordinates, 

= { Z1 B X + z 2 B_ 2 , (z x ,z 2 ) e C 2 } (4.19) 
and the action of SO (2) is then defined by 

(f>-(w,z) = (e^,e" 2 ^). 
With the notations of Section 4.3, we therefore have 

ipi = 4> and ci = —2. 
The expression for the cubic G-equivariant terms is 

C(M (3) ) = b l tr(MM)M + b 2 tr(M 2 )M + 63 (m 2 M + MM 2 - ^tr(M 2 M)Id \ ; 

with bj real coefficients depending on the specific Hamiltonian at hand. Setting M = Z\B\ + 
z 2 B- 2 we obtain after calculation in V^ 3 the expression 

C(M^) = 4 (V + |)|zx| 2 + {h + 6 3 )M 2 ) z 1 B 1 + 4 ((61 + b 3 )\ Zl \ 2 + h\z 2 \ 2 ) z 2 B_ 2 . 

Let us now check whether the 1x2 matrix A of Theorem 4.2 has maximal rank. From the 
above expression we deduce that 



en - c 2 i 



-26, 



C12 - c 2 2 



4h 



Therefore the maximality hypothesis is satisfied iff 63 ^ (which is a generic condition). 
Cases 4 and 5 are beyond the range of applicability of Theorem I 
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5 Appendix 



5.1 On the invariance properties of the resonance subspace 

In what follows we will sketch the proof of some of the facts about the invariance properties 
of the resonance subspace mentioned in the preliminaries section when (V, lo) is a symplectic 
representation space of the Lie group G and the Hamiltonian vector field A is G-equivariant. 

The resonance subspace U Ua is G— invariant Let A = A s + A n be the Jordan-Chevalley 
decomposition of A. Since by hypothesis A is G-equi variant, if $ : G x V — > V denotes the fr- 
action, for any g € G, we have that & g A = A& g . Equivalently, <& g A s + § g A n = A s & g + A n & g , 
and hence $ g A n & g -i + ® g A s & g -i = A n + A s . Since ® g A n <f> g -i is nilpotent, $ g A s $ g -i is 
semisimple, [$> g A n $ g -i , Q g A s <& g -i] = 0, and the Jordan-Chevalley decomposition is unique, 
we have that 

QgAn^g-! = A n and $ g A 3 $ g -i = A s , 

necessarily. This implies the G-invariance of U Vo = ker(e AsTv ° — I). Indeed, let v € U Vo . Hence, 
e AsTu °v = v. At the same time, for any g E G, 

e AaTvo (<& g v) = <& g e AsTvo v = $ g v, 

hence & g v G U Uo , that is, U Uo is G-invariant. 

The S 1 action and the G— action on U Uo commute Let : S 1 x U Vo — > U Vo be the 

S^-action on U Uo . For any g G G and any 8 E S 1 : 



<5> g ^e = $ g e eAs = e eAs <S> g = 



as required. 



5.2 Proof of Lemma 3.7 



The defining relation ( |3.4p of the function «i implies that for any vq E Vo, a, A G R, and £ G f) 
we have that 

(I - P) V^ o (fo - J 1+Q - )(«o + ui(«o, «, A, 0) = 0. 
Consequently, for any w% G V\\ 

(V^ - J 1+Q - K*)(fo* + ^(tuo, a, A, 0),m) 



d 
~d~t 

d 
dt 



t=o 

d(h x - 3 1+a - Kt)(tv + Vl (tv , a, A, £)) • Wl 



t=o 

d 2 (/^ - J 1+a - K*)(0)(u + £Vo«i(0 s a, A, • vo, w x ) 
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If we use ( |2.6| ) and ( |2.15 ), the previous expression can be matricially expressed as 



(0,1^; 



<r(\)hn T(\)hn + (V>(A) - (1 + a))hn ~ £ 

T(X)I 2n ~ (V»(A) " (1 + a))hn + £ PWhn 



v 

D Vo vi(0,a,X,£) ■ v 
= wj [r(A)I 2n - (V(A) - (1 + a))hn + £ + /o(A)£V «i(O s a, A, £)] «Q. 

Given that the previous equation is valid for no matter what vq G Vo and it?i G Vi, we can 
conclude that 

n rn x ^ ^(A) — (1 + a) T r(A) £ 
L>vb^i(0, a, A, 4) = — J 2n 7TTl2n 



p(\) ^ n P (xr n P (xy 

as required. ■ 
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